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1 Introduction 

In the modem world insurance companies operate in a financial environment. In particular, 
they invest their reserves into various assets and this may add more risk to the business. To 
our knowledge, the first model of an insurance company investing its capital into a risky 
asset has appeared in the short note llT2l where the author provided arguments showing 
that the asymptotic behavior of ruin probability is radically different from that in the clas¬ 
sical Lundberg-Cramer model. A rigorous analysis in da confirmed the conjecture. Using 
Kalashnikov’s estimates for linear finite difference equations with random coefficients, 1^ . 
it was shown that, independently of the safety loading, the ruin is imminent with unit prob¬ 
ability when the volatility a of the stock price is large with respect to the instantaneous rate 
of return a (namely, when 2a/cr^ < 1), and the rain probability is decreasing as a power 
function, when the volatility is small (namely, when 2a > 1). For the model with expo¬ 
nentially distributed claims the exact asymptotic was found. Threshold case 2o/cr^ = 1 was 
studied in the papers l29t . ll30l where it was shown, using techniques based on a renewal 
theorem, that the rain is imminent with unit probability. The setting of IT 3 I and (2^ is of the 
so-called non-life insurance: the company receives a flow of contributions and pay primes. 
The rain occurs when a new claim arrives and its value is too large to be covered by the re¬ 
serve: the risk process exits from the positive half-axis by a jump. The model can be studied 
in the discrete-time framework but the method of differential equations happens to be more 
efficient in the case of exponential claims where it allows to get the exact asymptotic. 

In the present note we investigate the setting of life annuity, or pension, insurance, when 
the company pays to the policyholder a rent and earns a premium when the insured person 
dies, see, e.g., 8 or l32l : in the classical literature such models are called models with 

negative risk sums. The highly stylized model of the collective risk theory may describe the 
situation where annuitant, entering the contract, pays a lump sum, e.g., the savings during 
his working life, getting in exchange regular payments until the death. One can interpret 
the premium as the reserve release. If the company has a large ’’stationary” portfolio of 
such contracts, with new arriving customers, then it is reasonable to think that the dates of 
reserve releases form a Poisson process and the sizes of premiums are independent random 
variables of the same distribution. A special type of life annuity, called le viager exists in 
France. This practice is, basically, an exchange of a property, a house or an apartment, for 
life-long regularly payments. The annuitant enjoys a regularly income without having to give 
a house. For the company having a portfolio of viagers the amount of premiums depends on 
the real estate market. 

For this classical model, where the reserve process has jumps upwards and may leave 
the positive half-axis only in a continuous way, the rain problem can be easily reduced to the 
rain problem for the non-life insurance using the so-called duality method, |^. Its idea is 
to define the ’’dual” process by replacing the line segments of the original process between 
consecutive jumps by downward jumps of the same depth and the upward jumps by line 
segments of the same height with positive fixed slope. Note that in the literature the models 
with upward jumps are often referred to as dual models, O, fll . 

The duality method does not work in our setting where the capital of the company or 
its fraction is invested into a risky asset. The change of two signs to the opposite ones 
in the equation defining the dynamics of the reserve leads to technical complications. In 
particular, now the rain may happen before the instant of the first jump and the latter is no 
more the instant of the regeneration, after which the process starts afresh provided being 
strictly positive. 
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Nevertheless, a suitable modification of arguments of fT3l and E9l combined with some 
ideas allows us to obtain the asymptotic of the ruin probability as the initial reserve tends to 
infinity. Expectedly, it is the same as in the non-life insurance case. In many countries there 
are rules allowing insurance companies to invest only a small share of their reserves into 
risky assets. Our simple model confirms that it is reasonable and even provides a quantitative 
answer. To avoid the situation when the ruin happens with probability one, the proportion of 
investment into the risky asset should be strictly less than 2a/(j^. 

It should be emphasized that the life insurance case is rather different and its study is not 
a straightforward exercise. The main difficulty in deriving the integro-differential equation 
is to prove the smoothness of the ruin probability and the integrability of derivatives. This 
issue is already delicate in the non-life insurance case. Unfortunately, it was not discussed 
it in GSi where the reader was directed towards the literature. Now we are not sure that we 
disposed at that time a clean reference. The smoothness of the exit probability is discussed 
in many papers, see, e.g., an interesting article f33t where it the explicit formula for an 
exponential functional of Brownian motion due to Marc Yor is used. Unfortunately, the 
needed smoothness property was established only under constraints on coefficients. One of 
the authors of the present note had a fruitful discussion with Marc on a possibility to deduce 
smoothness of the ruin probability function and the integrability of its derivatives without 
using complicated explicit formulae. His suggestions are realized here for the life insurance 
case. Note that in the literature on non-life insurance one can find other methods to establish 
the smoothness, for example an approach similar to verification theorems in the stochastic 
control theory, see (5). 

The structure of the paper is as follows. Section 2 contains the formulation of the main 
results. In Section 3 we establish an upper asymptotic bound for the exit probability (from 
]0, oo[) for a solution of a non-homogeneous linear stochastic equation and a lower asymp¬ 
totic bound for a small volatility case. As a tool we use the Dufresne theorem rather than 
Goldie’s renewal theorem from oa, which was the kew ingredient of arguments in ||2^ . 
The proof of Theorem |2] asserting that in the case of large volatility the ruin is imminent is 
given in Section 4. The regularity of the non-ruin probability ‘I’ is studied in Section 5 using 
a method based on integral representations. At the end of this section we derive the integro- 
differential equation for $. Section 6 contains the proof of the main theorem. Finally, in the 
appendix we provide a formulation of an ergodic theorem for an autoregression with random 
coefficients proven in f29l . 

Kalashnikov’s approach (developed further in his joint work with Ragnar Norberg (2T|) 
plays an important role in our study. Our techniques is elementary. More profound and 
general results can be found in ll27l . ll26l . f28l . ll23l . (24), ll25l et al. 

2 The model 

We are given a stochastic basis (17, F = (J^t)t>0i P) with a Wiener process W indepen¬ 
dent of the integer-valued random measure p{dt, dx) with the compensator p(dt, dx). 

Let us consider a process X — of the form 



(1) 


where a and a are arbitrary constants and c > 0. 

We shall assume that p{dt, dx) = adtF{dx) where F{dx) is a probability distribution 
on ]0, oo[. In this case the integral with respect to the jump measure is simply a compound 
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Poisson process. It can be written as where A'" is a Poisson process with intensity 

a and are random variables with common distribution F; the random variables W, N, ^i, 
i € N, are independent. We denote by Tn the successive jumps of N\ the inter-arrival time 
intervals — Ti_i are independent and exponentially distributed with the parameter a. 

In our main result (Theorem [T]l we assume that F is the exponential distribution with 
parameter 

Let ~ infjt : < 0} (the instant of ruin), <F{u) ~ P{t'^ < oo) (the ruin 

probability), and <P{u) := 1 — F{u) (the survival probability). 

The parameter values a = 0, cr = 0, correspond to the life insurance ( i.e. the dual) 
version of the Lundberg-Cramer model for which the risk process is usually written as 

Nt 

n :=u- + (2) 

i=l 

In the considered case the capital evolves due to continuously outgoing cash flow with rate 
c and incoming random payoffs at times forming an independent Poisson process N with 
intensity a. For the classical model with positive safety loading and F having a “non-heavy” 
tail, the Lundberg inequality provides an encouraging information: the ruin probability de¬ 
creases exponentially as the initial capital u tends to infinity. Moreover, for the exponentially 
distributed claims the ruin probability admits an explicit expression, see (2] or m. 

The more realistic case a > 0, a = 0, corresponding to non-risky investments, does not 
pose any problem (see, e.g., [m for estimates of the exit probabilities covering this case). 

We study here the case a > 0. Now the equation (|T]l describes the evolution of the 
reserve of an insurance company which pays a rent and continuously reinvests its capital 
into an asset with the price following a geometric Brownian motion. The same model can 
be used for the description of the capital of a venture company funding R&D and selling 
innovations, 141. 

Notations. Throughout the paper we shall use the following abbreviations: 

1 2 . 2k: 2a , ,,, 

K := a - -a , ^ - 1, rji ■-= K,t + aWt- 

2 

The solution of the linear stochastic equation (|T]| can he written using the Cauchy for¬ 
mula: 

Xt = e^*iu+ [ e~^^drs] . (3) 

V 4[o.t] J 

Theorem 1 Let F{x) = 1 — e~^l^, a; > 0. Assume that a > 0. 

(i) If P > 0, then for some K > 0 

Fiu) = Ku~^{I + o{l)), M —^ oo. (4) 

(ii) If P < 0, then 'F{u) = 1/or all u > 0. 

The formulation of this theorem is exactly the same as in fTSll for the non-life insurance 
model (the case P = 0 was analyzed in f29]| . (301 1. 

One must admit that in the life insurance models, especially, in the case of a company 
keeping a portfolio of viager contracts, the hypothesis that the benefits (i.e. prices of houses) 
follow the exponential distribution, are highly unrealistic. Nevertheless, we can claim, with¬ 
out any assumption on the distribution, that the ruin probabilities lay between two power 
functions, see Propositions [T]and|2l 
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The next result (implying the statement (ii) above) says that for 5 > 0 the ruin is immi¬ 
nent. It requires only that the distribution F has moments of positive order. 

Theorem 2 Assume that there is 8 > Q such that < oo. If /3 < 0, then Fiu) = 1 for 
any m > 0. 

The same model serves well in the situation where only a fixed part 7 g] 0, 1] of the 
capital is invested in the risky asset: one should only replace the parameters a and a in 0 
by ay and ay. Theorem [T] implies that the ruin with probability one will be avoided only if 
2ayl{ay'f' > 1 , i.e. when the share of investment into the risky asset is strictly less than 
2 a/(j 2 . 

It is worth mentioning that our conclusion are robust and holds for more general models. 
The reader may contest that intensity of outgoing payments c is constant. Indeed, after the 
death of the annuitant the payments stop and the intensity must decrease while with a new 
customer it increases. Easy comparison arguments show that the above statements hold in 
a more realistic situation where c = (cj) is a random process such that 0 < Ci < c < C 2 
where Ci and C 2 are constants. 

The crucial part of the asymptotic analysis in Theorem 1 is based on the fact that for 
the Markov process given by 0 the non-exit probability ${u) is smooth and satisfies the 
following equation: 

1 

—a‘^'L?I?''{u) + {au — c)^'{u) — a${u) + a I <l>{u + y)dF{y) — 0. (5) 

2 Jo 

With cr > 0 this equation is of the second order and, hence, requires two boundary 
conditions — in contrast to the classical case (a = 0 , a = 0 ) where it degenerates to an 
equation of the first order requiring a single boundary condition, see IIT 6 I . The estimate 
given in Proposition [T] shows that <P{oo) — 1 . 

There is an extensive literature concerning of the regularity of the survival probability 
for process with jumps in the context of non-life insurance models and models based on 
the Levy processes, see, e.g., (33), m, 03, ii, la. Our Theorem requiring only the 
smoothness of F and the integrability of F” seems to be the first result on the regularity of 
the survival probability in the considered setting. 


3 Asymptotic bounds for the small volatility case 

3.1 Upper asymptotic bound 

Let us consider the exit probability problem for a more general process X = of the form 

Xsds + a [ XsdWs-ct +Zt, (6) 

Jo 

where a,a,c > 0 are arbitrary constants and Z — {Zf)t>o is an increasing adapted cadlag 
process starting from zero. 

Proposition 1 Assume that in the general model 0 the parameter /3 > 0. Then 


Xt = u + a 


f 
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Proof. Let y be a solution of the linear stochastic equation 


Yt 


u + a 


[ Ysds + a [ 

Jo Jo 


Ys dWs - ct. 


Introducing the notation 


we express the solution as 


Rt := c f e dv 

Jo 

Yt ■.= e^^{u-Rt). 


(7) 


The difference X — Y satisfies the linear equation with zero initial condition. Since Z is 
increasing, we have the inequality X > Y showing that the exit of X from ]0, oo[ implies 
the exit of Y. Thus, 

Ip'(u) < P(i?oo > u). (8) 


and the asymptotic behavior of the ruin probability for this general model can be estimated 
by the tail behavior of the distribution function of Roo. Using the change of variable v = 
(4:la‘^)t and observing that Bt —(l/ 2 )crVL( 4 /^ 2 )( is a Wiener process we obtain the 

representation 


R 


OO 



^—{a—a‘^j2)v—(7Wv 


- i£ /'°° -2/3i+2Bt 

-Vo 



(-/ 3 ) 


The Dufresne theorem (see, O or d, Theorem 6.2) claims that is distributed as the 

random variable 1 /( 27 ) where 7 has the gamma distribution with parameter p. Thus, 


P{Roo > u) = P(2c/((t^7) > u)) = P (7 < 2c/((t^u)) = / dx 

1 [P) Jo 

2^ -/3 

~ a'^Pprip)'^ 


and the result follows from ([ 8 ll. □ 


3.2 Lower asymptotic bound 

The next result shows that the ruin probability decreases as the initial capital tends to infinity 
not faster than a certain power function. 

Proposition 2 Assume that P > 0. Then there exists /3* > 0 such that 

liminf T>{u) > 0 . 

U—¥00 


Proof. Let Y = (yfe)fe>i be the imbedded discrete-time process, that is the sequence of 
random variables defined recursively as follows: 

Yk = MkYk-i + Qk, Yo = u, (9) 


where 


rTk 

Mfc = and Qk = ~ c / dv . 

■ITk-i 
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Let 0“ := inf{fc : Vj, < 0}. It is clear that XT^. = Lfc for any k>l. So, for any u > 0 

= P(r“ < oo) > P(6l“ < oo). (10) 

Take £i g] 0, 1[ and chose B sufficiently large to ensure that 


Bi=B- 


> 0 . 


- q) 

Define the sets 

A = {Mk < £»} n {Qfc < e , Dfe = {Mfc < £1 n {Q^ < —B}. 

On the set Bf^ we have 

n n n 

Yn = uYlM, + J2Qk n + 


j—1 k—1 j—k-\-l 

Therefore, on the set B^, n D„+i 


Yn+l — Mn+lYn + Qn+1 < UQ + 


£<^(1 - q) 


£'(1 - q)' 


Ti BL — 1 j~) 

— B = ug — Bi. 


It is easy to check that itg” ^ < Bi, when u > Bi and 


n = 3 + 


ln(M/Bi) 


|ln£>| 


where [...] means the integer part. Therefore, 

P(r < oo) > P (n^=i A n D„+i) = (P(ri))" p(Di). 
Taking into account that P(A) > 0 and P(Di) > 0, we obtain that 

lim < oo) = oo 

ti—>-oo 

for any 


/3* > 


InP(ri) 

In g 


This implies the claim. □ 


( 11 ) 


4 Large volatility: proof of Theorem|2] 

We consider separately two cases: /3 < 0 and /3 — 0 and show that in both cases the ruin 
probability is equal to one. 

Proposition 3 Assume that /? < 0 and E < oo for some 5 > 0. Then B{u) = 1 for all 
u > 0. 
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Proof. As in the proof of Proposition |2] we consider the imbedded discrete-time process 
y = (Xk)k>i defined by (O- By virtue of ([TO} it is sufficient to show that P(0“ < oo) = 1. 
Note that for 5 g]0, —/3[ 


_ -£tg^S{Kt+<7Wt) _ ^St(^P+S)cT^/2 ^ 


and, therefore. 


EA/f = 


/*00 

/ ae~°‘*dt < 1. 

Jo 


( 12 ) 


According to (®,p. 250) if a random variable u is independent of W and has the exponential 
distribution with parameter a, then 


Ee 


q max^Tj<^iy{[i.v+Wv) _ /r 

-^20+7? — fi — q 


(13) 


provided that 


■\/2a + — q > 0. 

Changing the variable and estimating the integrand by its maximal value we get that 


E 


pOO 

Jo 

t>0 

_ C2^8{\n{26/a.) — l) ^ (/*u+crWi,) 

where v' is an exponential random variable with parameter a/2. In view of the equality ( 113b 
S (kv+(tWv) y/ (J^QL + K 

V (T^a + Ik? — Hi — Sa"^ 


provided that 


5 < 


y/a'^a + K? — K 


i.e. for such 5 > 0 we have the bound 

rTi 


Wo ^ 


'a + — K 


(T^a + — K — Sa'^ 


(14) 


Using these estimates and the assumption of the proposition we conclude that E IQi |'^ < oo 
for sufficiently small 5 > 0. Thus, the hypothesis of the ergodic theorem for autoregres¬ 
sion with random coefficients is fulfilled (see Proposition!^. Th^ latter claims that for any 
bounded uniformly continuous function / 


N 






( 15 ) 
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where 


OO k — 1 


c=Qi + '^2 n 

k=2 j=l 


Let us represent ( in the form 


C = ei- + 

Jo 


OO k—1 

Cl := E Qfc n 


fc=2 i=2 


Clearly, the random variables and (i?TiiE^ dv) are independent. Moreover, 

Lemma [T] given after the proof implies that the support of conditional distribution of the 
integral fg e^*~‘^'’dv given rjt = x is, unbounded from above. From this we easily infer that 
the support of distribution of is unbounded from below. Thus, for 


f{x) — l{a;<-l} + \^\ 1{-I<a;<0} 

the right-hand side of dll is strictly positive and, therefore, P(mffc>i Yj, < O) = 1. □ 
Lemma 1 Let ct > 0. Then the support of conditional distribution of the random variable 

T (tWv , 

7 = / e dv 

Jo 

given Wi = y is unbounded from above. 

Proof. It is well-known (see, e.g., 131]) that the conditional distribution of the Wiener process 
{Wv)v<i given Wi = y coincides with the distribution of the Brownian bridge with 
= Wt+t{y—Wi). Thus, the conditional distribution of I is the same as the unconditional 
distribution of 

Since the Wiener measure has the full support in the space Co[0,1] of continuous functions 
on [0,1] with zero initial value, the support of distribution of / is unbounded from above. □ 

Proposition 4 Assume that /? = 0 and E < oo for some <5 > 0. Then 'T{u) = 1 for all 
u > 0. 



Proof. In the considered case the imbedded discrete-time process is defined by ® with 

M/j = and Q/j = Cfe ~ c f g<r(WT^-Wv) 

Jtu-i 

where 14 = and AVf^ ;= I 4 —I4_i. To study the asymptotic properties of the equation 
(|9ll we use the approach proposed in (2^ for the non-life insurance models. To this end, 
define recursively the sequence of random variables putting 9q := 0 and 

Bn ~ mf{k > 6n-i: V/, - Vg^_^ < 0}, n > 1. (16) 

Note that = Y2i=i ‘^^3 where is a sequence of i.i.d. random variables dis¬ 

tributed as 9i. It is known (see, e.g., XIL7, Theorem la in El) that 

C := sup 

n>l 


P(0i > n) < OO . 


(17) 
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Putting t/fe = Yg^, 

we obtain from @ that 


Vk = akVk-l+bk, yo=u, 

where 

ASk 


flfc = n ^Sk-i+3 = 


i=i 

and 

A9k / A9k \ 


=E n ^6k-i+j Qok^ 


1=1 \j=i+i 1 


(18) 


It is clear that a^. < 1 a.s. Moreover, the first condition in Theorem|2]and the inequality ill4b 
with K = 0 implies that E|Qi 1^ < oo for any sufficiently small 5. Now, taking into account 
that 


Ae, 


ABi / ABi \ 

|6iI<E n IOll = E 

1=1 \j=l+l I 


ai 


A01 


\Qi\ < El'SiN 


1=1 11^=1 1=1 
we can get, for r £]0,1[ and an increasing sequence of integers In, that 

E|6ir<l + r^ 1 P(|6^|>^) 


n>l 


n>l \J —1 / n>l 

<l + rE|Qi|^^-^+rC^^ 


n>l 


n>l 


1/2 ■ 


Putting here In = [n^^] we obtain that E < oo for any r g] 0, 5/5[. Therefore, due 
to Proposition [8] we obtain that for any bounded uniformly continuous function / 


N 




fc=l 


oo k — 1 


(19) 


where 

c=6i+E n 

k=2 j=l 

Now we show that 

P(C < —x) > 0 for any a; > 0. 

Indeed, the random variable ( I19l l can be represented as 

OO k—1 

c = fei+oiCi, Ci = E^fcn“i' 

fc=2 j=2 

It is clear that the (^i is independent on bi and ai. Note that on the set {A9i = 1} we 
have fli = Ml and bi = Qi - Therefore, for any x > 0 

P(C < —x) > P(6i + ai (^1 < —X, Adi = 1) = P(Qi + Mi( < —X , wt^ < 0) 


and we conclude as in the previous proposition. □ 
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5 Regularity of the ruin probability 

5.1 Integral representations 

The proof of smoothness of a function H admitting an integral representation is based on a 
simple idea which merits to be explained. 

First, we recall the classical result on differentiability of the integral Hiu) = f f(u, z)dz 
where f{u, .) € for each u from open subset [7 G M. If z) is differentiable on an open 
interval ]mo — e, uq + ^[for almost all z and, on this interval, \df{., z)/du\ < g{z) (a.e.) 
where g G L^. Then H is differentiable at mq and H'[uf) = f df(uQ, z^jdudz. 

Suppose that we are given a bounded measurable function h(z) and a Gaussian random 
variable C, ~ Af(0,1). Let H(u) = Eh(u + 0 = J h{u + x)<pop{x)dx. Then H is differen¬ 
tiable and even of the class (7°°. Of course, the above result cannot be applied directly. But 
using the change of variable we get the representation 

H{u)^ I h{u + x)tpQ^i{x)dx = j h{z)(pQp{z — u)dz. 

Now the parameter u appears only in the function the integrand is differentiable in u 
and we can apply the classical sufficient condition. 

The issues here are: an integral representation, the smoothness of the density, and the in- 
tegrability of its derivatives. In the case of survival probability <P the integral representation 
is obtained from the strong Markov property. Unfortunately, the structure of the represen¬ 
tation is rather complicated, the random variable standing for ^ is not a Gaussian one, and 
its density is not given by a simple formula. Nevertheless, the idea of using a change of 
variable to move the parameter from the unknown function on which we have only a limited 
information (essentially, the boundedness and measurability) to the density still does work. 
The main difficulty is to check the smoothness of the density and find appropriate bounds 
for its derivatives. 


Theorem 3 Assume that the distribution function of^i has a density f differentiable on M-|- 
and such that f' G Z/^(K+). Then T>(u) is two times continuously differentiable on ]0, oo[. 

Proof. We again consider the process 


= e"* (u - Rt), (20) 

where Rt is defined in ([7]l. Put 

0“ := inf{f >0: U“ < 0}. (21) 

By virtue of the strong Markov property of X“ 

<?(u) =E<I>(Xe“»ATj- (22) 

Note that the process y“ is strictly positive before the time zero at 0“, and strictly 
negative afterwards. Due to the independence of the Wiener process and the instants of 
jumps, 0“ Ti a.s. Thus, > 0} = {0“ > Ti} a.s. Taking into account that = 0, 
we get that 

>Piu) = E/|y„ >0}^(2fTj = +Cl) = <Pliu)+4i2iu), 
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where 


^i(m) 



^2(w) := ot 



EG(y“)e““*dt 


with 

G{y) ■- l{y>o}E<?(t/ + Ci) = 

We analyze separately the smoothness of <l>i 
integral representations. 


pOO 

^{y>o} / Hy + x)dF{x). 

Jo 

and <p 2 using for this function appropriate 


5.2 Smoothness of $2 


We start with a simpler case of $2 and show that this function is infinitely differentiable 
without any assumptions on the distributions of ^ 1 . 

From the representation 




-j: 


^ds, 


we obtain, using the independence of y“ and the process ( 77 s — Vi)s>i^ that 

E(G(yt“)|yi“) = G(f,Fi“), 

where 


G{t,y)--EGl^e^^-^^y-c j e^*~'^‘dsj. 

Substituting the expression for given by ( 122b we have: 

<? 2 (m) = ® ^ E(G(y“)|yr)Q;e““‘df = EH (u - i?i)) 

where 7F is a function taking values in [ 0 , 1 ] and given by the formula 


H{y)--a j G(t,y)e °‘*dt. 

Taking into account that the conditional distribution of the process (iys)s<i given lEi = x 
is the same as of the Brownian bridge = {Bs)s<i with = Ws + s{x — Wi) we 
obtain the representation 


‘^ 2 ( m ) = I EH{e'^~^'^'^{u-0)ipo,i{x)dx 


where 


C := c 


f 


-{ks-\-sx-\-(t{Ws-sWi)) 


(23) 


(24) 


Lemma |2] below asserts that for every x the random variable has a density p{x ,.) on 
]0, oo[ and we easily obtain from ( 123b by changing variable that 

<^ 2 (u) = J J H{e'^^'^^z)p{x,u — z)ipQp{x)dxdz. (25) 
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Lemma 2 The random variable C,^ has a density p{x ,.) € C°° such that for any n > 1 

sup j—p{x,y) ( 26 ) 

y>0 <JV 

with some constant Cn and {cf^ /dy^)p{x, 0) = 0. 

Proof. We obtain the result using again the integral representation. Let us introduce the 
random process 

Ds := [{Ws — 2 sWi/2) + s{Wi/2 — W^l))-^{s<l/2} 

+ ((1 - s)(W. - W,/2)) - s(Wi - Ws))^s>l/2}, 

and the piecewise linear function 

7s := •s/{s<i/2} + (1 — s)-^{s>l/2}i s ^ [0, !]• 

The following identity is obvious: 
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and, similarly, 

an 

— Q^°\x,vix,y)) 

It is easily seen that 


^\x,v{x,y)) 

K{x,v{x,y)) 


Q^'^\x,v{x,y)). 




V’0,l/2(z} 
it'"+i(a;, z) 


^ ^ (^)^n—A: (^5 i 

k=0 


(28) 


where Pk{z) is a polynomial of order k and Rn_f^{x, z) is a linear combination of products 
of derivatives of K{x, z) in variable z. Note that for any x,y £ [0,1] we have the bounds 


— \x\ — (j\z\ — ScrW^i < Ks + sa; + aDs + o'ysZ < k + |a;| + cr| 2 | + SaW* 


where = supg <2 | Ws|- It follows that there exists a constant Cn > 0 such that 

\Q^'^Hx,z)\ < + (29) 

Since E e^” < oo, for each x, y and n the derivative {d^/dy^) {x,v{x, y)) admits 

a P-integrable bound. Thus, we can differentiate under the sign of expectation and obtain 
that 


^ p(x y) — E ^ ^ d q( 0)/ / NN EQ^"Va; v(x y)) 

Qy„P[x,y) i^Qyn K{x,v{x,y)) dy^^ [x,v[x,y)) [x,v[x,y)). 


Moreover, the bound ll29b ensures that, for some constant Cn 


supE 

v>o 


an 

— Q^°\x,vix,y)) 


< E sup I(a;, a)I < 
zeR 


and the bound ( 126b holds. 

Since v{x, 0+) = oo, the bound ( 129b implies that /dy'^)p{x, 0) = 0. □ 


Remark 1 It is worth to trace in these arguments the dependence of the constant 67" on the 
parameters c and a when they are approaching zero. From the formula ( 128b it is clear that 
Cn should be proportional to 


Proposition 5 The function <p 2 {u) belongs to 67“ (]0, oo[). 


Proof. Putting 

H{u,z) ■.= j z)p{x,u — z)tpop{x) dx, 

we rewrite the formula (125b as 


L 


^ 2 («) = / H{u,z)dz. 


( 30 ) 


Clearly, the function H is continuous on ]0, oo[x]0, oo[. Using Lemma|2]we obtain that 


du 


!■ 


d 


H{u,z) = I H{e ^ z)—p{x,u-z)ipop{x)dx 
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and 


sup —H{u,z) < oo. 

U,Z 


By induction, for every n > 1, 



and 


9 ” ~ 

U,Z 


By virtue of Lemma|2p(a;, 0) = 0, i.e. H{u, u) = 0. So, 



In the same way we check that 



for any n > 1. □ 

5.3 Smoothness of $i 

Arguing in the same spirit as in the previous subsection but taking this time the conditional 
expectation with respect to Wt we obtain that 



where we use the abbreviations 


h{u,t,x) := h{y) = E${y + ^i) 


and 



It is easily seen that the random variable has infinitely differentiable density (the 
same as of defined in (I24ll but with the parameters ct, Kt, and Unfortunately, 

derivatives of this density have non-integrable singularities as t tends to zero (see Remark at 
the end of previous subsection). By this reason we cannot use the strategy of proof used for 
4>2. Nevertheless, the hypothesis on the distribution of allows us to establish the claimed 
result. 

Note that the function x —>■ is strictly decreasing and maps K onto M+. Let denote 

z{t ,.) its inverse. The derivative of the latter is given by the formula 
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where 


L{t, z)=c [ se-^^z/t+i-s+aiW.-{s/t)W,) 

Jo 

Changing the variable we obtain that 

I ^{C’^<u}Hu,t,z)ipo,i J Ku,t,z{t,x)) D{t,z{t,x))dx, 

where 

Summarizing, we get that 


(31) 


L{t,z) 


$i{u) = aE j 
Jo 


— at 1 , 

e at, 


where 


H{t,u) := [ h(u,t, z(t,x)) D(t, z{t,x)} dx. (32) 

JO 


Proposition 6 Under the conditions of Theorem\^the function H{t,u) defined in mfor 
any fixed uq > 0 satisfies the following inequality 


sup E sup (|//u(t,u)| + m)|) < oo. 

tG]0,2] u>Uo 


(33) 


Proof. In virtue of the hypothesis the function h{y) = E<P{y + ^i) is differentiable. Differ¬ 
entiating ll32t we get that 

Hu{t,u) = h{0) D{t, z{t,u)) + [ hu{u, z{t,x)) D{t, z{t,x)) dx, (34) 

Jo 

where h(0) = E<l>(^i) and 

hu{u, t, z{t, x)) = h'[u - x)). 


Note that 


d 


— h(u,t,z(t,x)) = hu(u,t,z(t,x))[azx(t,x)(u — x) — l]. 


Therefore, 

hu{u, z{t, x))D{t, z(t, x)) dx = — 


pv 

Jo 


pu 

Jo 


ipo,l {z{t,x)/Vi) 


Q at{u —x) + L{t,z{t,x)) 
Integrating by parts and taking into account that z{t, 0-I-) = oo we get that 


dxh(u,t,z{t,x)). 


H. 

where 

0[u, t, z) = 


ru 

(t,u) = / h{u,t,z{t,x))0{u,t,z{t,x))(po^i (^z{t,x)/y/i'j dx, (35) 

Jo 


taL{t, z) + tLz{t, z) 


L{t, z){at {u - C’^)) + L{t, z)) Li^t, z) (at {u - C*’^) -b L{t, z)) 


2 • 
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Inspecting the formula ( I31t defining L{t, z) we conclude that there exist positive con¬ 
stants Co (’’small”) and Ci (’’large”) such that 

max {L{t, z) + t\Lz{t, z)|) < ^ 

te]o, 2 ] 

and for any t g]0 , 2] 

where Wt ~ max^,<t |H4i|- Taking this into account we obtain that for some C > 0 

\0{u,t,z)\ < . (36) 

Using a generic notation C for a constant (which may vary even within a single formula) we 
obtain, for any t €]0, 2] 

pCO 2 

\Hu{t,u)\ / e^^"^~^L{t,z)dz 

J z{t,u) 




f 




So, we have the bound 

For any u > uq > 0 and t €]0, 2] 




^2 

T dz. 




i.e. 

1 

\z{t,u)\ > - In --- w*. 

a tc 

Put 

to := 2^ , F ■.= {Wt* < 

Thus, for t g] 0, fo] on the set F we have the inequality 

\zit,u)\ > 1. 

Taking into account that E e“^* < oo for any a and I > 0, we obtain that 

E max sup \F[u{t,u)\ fCCfT^Ee^'^^ < oo. 
t^[to,2] u>0 

For t e]0, to] we have 

E max \Hu{t,u)\ < Ct~'^ fe"® -f Ee'^'^* IpA 
U^Uq \ J 


< Ct~ 


It —— 

' I ^ A4- 


+ 


VEe'^ty* y'p {w^ > ii/4) 
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By the Chebyshev inequality we have: 


P (Wt <e~^ =e“* 


that is 

supe* '' P (Wt > < oo. 

t>0 ' 

This implies that 

max E max \Hu{t,u)\ < oo. 
tS]0,to] u>uo 

Therefore, 


max E max \Hu{t,u)\ < oo. 
tG]0,2] u>uo 


Differentiating ( 135b we find that 


u) 


h(0) 0{u, t, z{t, u)) If 


z{t, u) 


)^i: 


r{u,t,z{t,x))ifi 


( z{t,x) \ 

\ Vi ) 


dx, 


where 

T{u,t,z) = hu{u,t,z)0{u,t,z) + h{u,t,z)0u{u,t,z). 

By assumption, the distribution function F has the density / whose derivative f' is a con¬ 
tinuous function on K+ integrable with respect to the Lebesgue measure. Changing variable 
we get that 


J rco 1 POO 

V{y) = ^ ^{y + x)f{x)dx=— F{x)f{x-y)dx 

dy Jo dy Jy 

poo 

= -^(y)/(0)- / ^{z)f'{z-y)dz, 


i.e. 

sup \h'{x)\ < oo . 

a ::>0 

Similarly to ( 136b we obtain that 


\Ttiu,z)\ < 


Therefore, 

\Huu{t,u)\ < e“^^, 

implying that 

max E max \F[uu{t,u)\ < oo. 
t6]0,2] u>uo 


Proposition[^is proven. □ 
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5.4 Integro-differential equation for the survival probability 


Proposition 7 Suppose that $ G C^. Then satisfies the equation (|j]). 


Proof. For h > 0 and e > 0 assumed to be small enough to ensure that u G]e, e ^ [ we put 

inf {t > 0: ^ [e , e“^]} A h. 

Using the Ito formula and taking into account that on the interval on [0, Ti [ the process X“ 
coincides with we obtain the representation 

ATi ) = M + n / <l>'{Yfi) dWs 

+ J + (ay“ - cy(y“))ds 

+ 6) - WjUm<r^}- 

Due to the strong Markov property 4>{u) = For every e > 0 the integrands 

above are bounded by constants and, hence, the expectation of the stochastic integral is zero. 
Moreover, A Ti — h when h is sufficiently small (the threshold below which we have this 
equality, of course, depends on cj). 

It follows that, independently of e. 


1 /I \ 1 

-E J [^rj^{Ys" f<P"{Y,^) + aY^$'{Y^) - c^s ^ -o^r?T>''{u) + (an - c)T>'(u). 


Finally, 


iE(<?(y^^+ei)-= [ \T>{Yt^ + ^i) - ${Yn)e-‘^Ut. 

•Z 0 

The right-hand converges to a(E^(w + ^i) — 'T{u)) as /i —>• 0 in virtue of the Lebesgue 
theorem on dominated convergence. It follows that 4> satisfies the equation (|5]l. □ 


6 Proof of Theorem [T] 

Assume that the claims are exponentially distributed, i.e. F{x) = 1 — ^. Similarly to 

the classical case, this assumption allows us to obtain for the ruin probability an ordinary 
differential equation (but of a higher order). Indeed, now the integro-differential equation 
© has the form 

1 

-a^v^$"{u)-\-{au — c)^'{u) — a.${u)-\-— / T>(u^dy — fS. (37) 

^ P Jo 

Notice that 

I pCO ■, pCO 

— I T>{u + y)e~^^^dy = —<T(u) H-/ $(u + y)e~^^^dy. 

du Jo P Jo 

Differentiating 07l i and adding to it the obtained identity multiplied by p, we exclude the 
integral term and arrive at a third order ordinary differential equation. It does not contain 
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the function itself and, therefore, is reduced to a second order differential equation for the 
function G = <P' which can be easily transformed to the form 


where 


G" - p{u)G' + po{u)G ^ 0, 

p{u) :=- 2 (l -\ - 2 )-f 

a \ (7^) u 


(38) 


a \ 1 2c 1 


, , 2a 1 , / N 2 1 

PoW :=- 2~ + («-« + c/p) -o 

u ( 7 ^ u 


The substitution G{u) = R(u)Z{u) with 

R{u) ■.= exp 




eliminates the first derivative and leads to the equation 

Z” - q{u)Z = Q, 

where 


(39) 


i=2 


dx < 00 , 


with certain constants Ai which are of no importance in our asymptotic analysis. It is easy 
to check that 

where xq = sup{a; > 1: q{x) < 0} + 1. 

According to fTOll . pp. 54-55, the equation ( 139b has two fundamental solutions and 
Z- 


f 

•J Xr\ 


Z±{x) = exp {*! y/q{z)dz^ (l-l-o(l)), 


00, 


Z± ~ exp 


Since 


R{x) ~ exp 


{§ 


a — a , 

-I-T— In® 


a -\- (T 


ln®| 


we obtain that ( 138b admits, as solutions, functions with the following asymptotics: 

G+(®) ~ G_(®) ~ . 

The differential equation of the third order for <P has three solutions: <Po{u) — 1, 


pu p-roo 

$+{u) = / G+{x)dx $-{u) = — / G-{x)dx. 

J Xq J U 

The ruin probability (Z/ := 1 — is the linear combination of these functions, i.e. 


>P(u) = Co + Ci<P+{u) + C2<P-iu). 
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Since <P+{u) oo as u —>■ oo, we obtain immediately that Ci = 0. For the case /3 > 0 we 
know from Proposition [T]that <f'(oo) = 0. Thus, Cq = 0 and 

r°° 2 

<P(u) = C 2 / (1 + 5{x)) dx, 

J U 

where 5{x) —>• 0 as a; —>• 00 and C 2 > 0. The integral decreases at infinity as the power 
function u~^ j/S and we obtain Theorem [T] □ 


7 Appendix: ergodic theorem for an autoregression with random coefficients 

Let {an ,bn)n>i be an i.i.d. sequence of random variable in and let xq be an arbitrary 
constant. Define the sequence (xn) recursively by the formula 

x„ = a„Xn-i + b„, n> 1. 

Proposition 8 Assume that there exists S (E]0, 1] such that 

E|ai|‘^<l, E|6i|^<oo. 

Then for any bounded uniformly continuous function f 

1 ^ 

n—1 


where 


00 k— 1 


^=z] n 

fc=i j=i 


with 


0 

n ^ ^ 


i=i 


The proof is given in f29ll . 
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